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LARGE  AVIATION  FUEL  FIRES 


EXECUTIVE  SUMMARY 


Evaluation  and  quantification  of  the  response  time  of  various  munitions 
to  accidental  fire  impingement  is  an  integral  part  of  the  DOD's  ordnance 
thermal  protection  program.  The  thermal  protection  program  involves 
interaction  of  numerous  thermodynamic  and  heat  transfer  disciplines.  This 
report  preserits  exclusively  the  efforts  on  numerical  and  experimental 
simulations  of  the  effects  on  ordnance  cookoff  time  when  external  cooling  is 
applied . 

To  assess  effects  of  the  external  cooling  on  fire-exposed  ordnance 
cookoff  time,  a  validated  heat  transfer  computer  model  was  adopted  and 
modified  to  include  principles  describing  the  external  cooling  effects.  The 
new  model  predicts  the  change  in  ordnance  cookoff  time  as  a  function  of 
coolant  application  rate.  Calculations  demonstrate  that  external 
cooling  can  delay  the  cookoff  time. 

The  thermal  interaction  between  the  incident  flame  and  ordnance,  as 
well  as  the  effects  of  coolant  on  cookoff  time,  were  simulated 
experimentally.  Specially  designed  calorimeters  were  instrumented  and 
placed  in  a  pool  fire  to  measure  the  transient  heat  flux  and  to  quantify  the 
effects  of  various  coolants.  Effect  of  an  Intumescent  coating  used  on  various 
ordnance  was  evaluated  experimentally  and  compared  with  the  response  of  the 
thermally  unprotected  ordnance.  The  limited  experimental  data,  although 
inconclusive,  suggest  that  thermally  unprotected  ordnance  can  achieve  longer 
cookoff  times  than  coated  ordnance,  when  external  cooling  is  applied. 
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SECTION  I 


INTRODUCTION 


A.  OBJECTIVE 

This  project  was  conducted  to  calculate  the  response  time  for  fire- 
exposed  ordnance  and  determine  how  long  ordnance  cookoff  can  be  delayed  by 
application  of  external  cooling.  Different  coolant  types  and  application 
techniques  were  tested  to  develop  pragmatic  criteria  for  response  times  and 
cooling  procedure-  to  be  used  for  ordnance  exposed  to  fire. 


B.  BACKGROUND 

The  US  armed  forces  have  had  a  number  of  incidents  in  which  ordnance 
exposed  to  file  tidh  cooked  off>  causing  majol  loss  of  life  and  property. 
Incidents  on  the  USS  Forrestal,  USS  Enterprise,  and  the  USS  Nimitz  al... 
resulted  in  losses  of  life  and  materiel  (References  1,  2,  and  3).  The  USaF 
has  also  had  such  incidents;  for  example,  a  maintenance  accident  at  Mountain 
Home  AFB  resulted  in  the  burning  of  an  F-111  loaded  with  Mark  82  bombs.  In 
this  incident  the  firefighters  retreated  and  all  of  the  bombs  cooked  off, 
resulting  in  significant  materiel  losses. 

In  the  event  of  a  postcrash  fire,  it  is  likely  that  ordnance  will  be 
exposed  to  high  heat  fluxes.  Ordnance  reactions  (propulsion  or  detonation 
depending  on  the  nature  of  the  weapon)  to  a  fire  environment  can, range  from 
mild  burning  to  violent  explosion.  The  extent  of  the  reaction  Is  determined 
by  Che  Intensity  and  duration  of  heating,  and  the  thermal  protection  of  the 
ordnance.  The  ordnance  may  have  a  cookoff  reaction  after  removal  of  the 
heat  source  (usually  by  fire  extinguishment)  because  of  self-beating  from 
internal  exothermal  reactions.  Firefighting  efforts  and  the  associated 
hazards  are  complicated  when  fire -exposed  ordnance  ere  present.  Criteria 
are  not  available  to  accurately  define  the  cooling/handling  requirements  and 
safe  response  procedures  to  be  used  for  ordnance  exposed  to  fire. 
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C.  SCOPE 

The  scope  of  this  project  was  to  Investigate  and  evaluate  all  existing 
ordnance  response  models  and  to  modify  a  viable  candidate  to  include  the 
effects  of  external  cooling.  The  modified  model  was  to  be  capable  of 
predicting  the  cooling  required  to  prevent  a  runaway  reaction  and  to  assess 
the  cookoff  time  extension  resulting  from  external  cooling.  Small-scale 
calorimeter  tests  were  conducted  to  obtain  essential  parameters  needed  to 
calculate  cookoff  time  extensions  which  can  be  achieved  by  application  of 
external  cooling.  To  evaluate  the  effects  of  intumescent  thermal  protection 
on  cookoff  time  extension,  coated  calorimeters  were  tested  and  the  data  were 
compared  with  the  response  data  from  uncoated  calorimeters. 

D.  APPROACH 


This  report  documents  numerical  modeling  conducted  as  one  component  of 
a  long  range  effort  to  reduce  losses  from  ordnance  cookoff.  Numerical 
models  which  predict  the  transient  response  of  ordnance  exposed  to  heat 
fluxes  have  been  developed  in  the  past.  One  of  these  models  has  been 
updated  to  include  external  cooling  and  its  effect  on  ordnance  response. 
Testing  of  subscale  models  has  been  used  to  measure  the  effectiveness  of 
various  coolants,  with  their  heat  transfer  rates  used  as  input  to  the 
computer  model.  The  change  in  cookoff  time  as  a  function  of  coolant 
application  rate  is  then  predicted. 
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SECTION  II 


MATHEMATICAL  MODELING 

A ,  INTRODUCTION 

The  derivation  of  heat  transfer  equations  describing  the  thermal 
loading  of  an  ordnance  is  based  on  the  physical  model  shown  in  Figure  1. 

In  this  figure  various  components  of  a  typical  ordnance  (a  Mark  82  general 
purpose  bomb)  are  shown  in  a  sector  cross  section.  Each  layer  is  characterized 
by  its  thickness  and  thermal  properties;  perfect  thermal  contact  between  each 
of  the  layers  is  assumed.  Internal  energy  generation  due  to  chemical  reaction 
is  allowed  in  any  layer.  The  energy  generation  is  modeled  by  the  use  of 
first-order  Arrhenius  kinetics.  The  ordnance  is  assumed  to  be  Immersed  in 
a  large  aviation  fuel  fire,  which  is  characterized  by  a  flame  temperature.  The 


Figure  1.  Sector  of  Cylindrical  Cross  Section  of  a  General-Purpose  Bomb, 
Mark  82  (extracted  from  Reference  4). 
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heat  flux  input  from  the  fire  to  the  ordnance  is  modeled  with  radiative  and 
convective  components.  Ordnance  cookoff  Is  defined  as  the  point  where  the 
explosive  begins  a  thermally  induced,  exothermic,  runaway  decomposition. 
This  process  is  modeled  by  use  of  the  first-order  Arrhenius  kinetics. 

B.  GOVERNING  EQUATIONS 

The  governing  equations  describing  the  overall  heat  transfer  processes 
in  a  fire-exposed  ordnance  must  include; 

1.  Equations  describing  the  incident  external  heat  flux. 

2.  Equations  describing  the  heat  transfer  through  various 
Internal  layers. 

C.  EXTERNAL  HEATING 


The  total  incident  heat  flux  on  the  exterior  boundary  of  the  ordnance, 
q^,  is  equal  to  the  sum  of  radiative  and  convective  heat  flux. 

Of  -  <  »  -  t‘)  *  h  <Tj  .  I^) 

with  the  terms  defined  as  follows: 


T^:  flame  temperature 

T^ :  surface  temperature 

t  ■.  emissivlty  factor 

a‘.  Stefan  Boltzman  constant 

h:  convective  heat  transfer  c'.eff Icient 

The  incident  flame  temperature  history  used  in  the  program  was  taken  from 
JP-5  fire  data  presented  by  Russell  and  Canfield  (Reference  4).  The 
emissivity  factor  of  the  fire  is  assumed  to  be  0.99.  The  convective  heat 
transfer  r oefflcient  is  taken  from  Boyer  and  Russell  (Reference  5). 

During  cooling,  the  surface  boundary  condition  is  modeled  using 
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q  -  U  (T  -  T  ) 
oa  w  s 

where  U  is  the  overall  heat  transfer  coefficient  and  T  Is  the  coolant 
oa  w 

temperature.  The  values  for  were  obtained  from  experimental  data 
provided  by  the  subscale  calorimeter  tests. 

D.  INTERNAL  HEATING 


Heat  is  transferred  through  various  layers  of  munitions  by  conduction. 
The  governing  internal  heat  transfer  equation  is  the  algebraic  sum  of 
conduction  and  internal  energy  of  the  materials  at  each  layer.  The  Internal 
heat  transfer  equation  in  cylindrical  coordinates  becomes: 


dTj 

‘"pi  JT 


-  k, 


& 

3r 


(  r- 


dr 


)  1  +  q. 


i  -  1 . n 


(1) 


The  terms  in  Che  energy  equation  are  defined  as  follows: 


i:  layer  index 

n:  total  number  of  layers 

3 

p:  density  (g/cm  ) 

C^:  specific  heat  (cal/g  K) 

T:  temperature  (K) 

t:  time  (seconds) 

2 

k:  thermal  conductivity  (cal/cm  s  K) 

r:  radius  (centimeters) 

q:  rate  of  Internal  energy  generation  (cal/s) 


The  rate  of  internal  energy  generation  (q)  is  expressed  as  a  first-order 
Arrhenius  reaction: 


q 


p-i  % 


z 

p 


exp 


-E  /  R  T  ) 
P  i 


(2) 


where 
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m:  number  of  components  of  the  explosive 
i:  represents  the  calculation  node 
Q:  heat  of  reaction  (cal/g) 

Z;  collision  number  (liters/s) 

E;  activation  energy  (cal/mole) 

R:  universal  gas  constant  (cal/mole  K) 

E.  BOUNDARY  CONDITION  AND  ASSUMPTIONS 


The  above  governing  equations  are  solved  according  to  the  following 
boundary  conditions: 


The  outer  boundary  condition  is 


k  4 

?r  -  * "  <  Tf  *  “"s  >  ^  ^  > 


(3) 


Similarly,  the  boundary  condition  for  the  period  of  cooling  becomes 


iT  “  "oa<  \  ^  > 


(A) 


The  boundary  condition  at  each  interface  assumes  that  the  heat  flux  across 
the  boundary  is  equal  on  both  sides  of  the  interface  and  the  temperatures 
are  the  same.  This  is  expressed  as: 


<3T,  , 

k  — ~  —  k  — 
i  dx  *^1+1  dr 


(5) 


and 


T 


i 


^e+1 


(6) 
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Furthermore,  it  is  assumed  that  the  variation  of  the  temperature  at  the 
center  of  the  ordnance  is  finite,  therefore 


ar* 


-  0,  at  r-0 


and  the  ordnance  is  initially  at  the  ambient  temperature 


T  =  T 


amb 


Y.  FINITE  DIFFERENCING  METHOD 


(7) 


(8) 


The  governing  heat  transfer  principles  and  the  boundary  conditions 
described  In  Equations  <1)  through  (8)  are  solved  numerically  using  a  finite 
differencing  scheme.  The  algorithm  used  to  obtain  an  approximate  numerical 
solution  is  based  on  an  explicit  central  differencing  method  (Reference  6). 
The  spatial  location  of  each  finite  element  is  based  on  the  cylindrical 
distribution  of  Ar  as  shovm  in  Figure  2. 


Expressed  in  the  f Inite -difference  form,  the  partial  differencial 
Equation  (1)  describing  the  energy  transfer  at  each  layer  becomes 

^  0.5  ■  ij ) 

i  i  r^ 

[Af|  ^  0,25  -  0.5  Ar^)[Tj  -  t|^^  ] 


Af|  (Ar)^ 


(9) 
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where  1  represents  the  current  node  and  j  represents  the  current  time.  The 
internal  energy  term  is  expressed  as 


where  n  represents  the  current  reactive  component,  m  represents  the  total 
number  of  reactants,  and  £  represents  the  current  layer.  The  Fourier 
modulus,  AF,  is  defined  as 


At 


(11) 


Figure  2.  Elemental  Surface  Layer  Analog. 
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The  surface  boundary  condition,  Equation  (3),  in  finite  difference  form  is 


^ 

‘  ^  1 

1 


q  +  2 


APJ  (r^)(Ar^) 

(r^  -  Arj^/4) 


[  h  (T^  - 


tJ)  + 


<7  (T^  -  (t|)^)]  -  Af| 


(r  -  Ar  /2) 

- - - i -  f  q'J  .  X- 

(r^  -  Ar,/4)  ^  "2 


(12) 


The  surface  boundary  condition  during  cooling  is  expressed  as 


^ 

"  "  ^  kJ 

1 


q  +  2 


AFj  (r^.)(Ar^) 

k|  (r^ 


Ar^/4) 


(  U  (T 
oa  w 


T|)  ]  - 


(r,_  - 

“'l  (t,  -  4t,/4)  ‘  H 


xJ  1 
2  - 


(13) 


The  interface  boundary  condition  between  layers  is  expressed  as 


xJ+1  _  tJ  +  _ 

1  1  Ax 

‘^l^  *  — ) 


Iv  h.0<^0  * 


28^  AXj^ 

Ax. 

K,(r  — ) 


Ax  T^  -  T^  28  AX 

•  •  2  ^  ~  - 


Ax, 


Ax 


Ki("-  — ) 


Ax  T^  -  T-l  28.  AXt 


Ax, 


Ax.  Ax. 

[(^1  +  4  )(1-)  <1 


Ax. 


(14) 


where : 


K^At 


1  /A  ^2 

p^c^(AXj^) 
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^1+1/2 


Ax 


\/2  ~  2 


1-0 


[  “  +  (  *  )  -  1)1 


The  coefficient  ^  is  the  gray-body  shape  factor  resulting  from  the 
radiative  heat  transfer  processes  (Reference  7).  Tlie  center  boundary 
condition  expressed  in  finite  difference  becomes 


- 


t|  +  2  ( 


k|  At 


Cp^  (At^) 


2  >  <  ■^1-1 


-  )  + 


Pi  Cp. 


(15) 


A  FORTRAN  program  was  written  to  solve  the  governing  differencing 
Equations  (9)  through  (15).  Details  of  the  computer  model  are  discussed  in 
the  next  section. 
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SECTION  III 


COMPUTER  MODEL 


A.  INTRODUCTION 

The  finite  difference  equations  used  to  calculate  the  temperature 
profiles  through  the  ordnance  were  implemented  in  a  computer  program.  A 
flowchart  of  the  program  is  contained  in  Appendix  A  and  the  program  listing 
is  in  Appendix  B.  The  program  is  a  modified  version  of  one  originally 
written  by  Russell  and  Canfield  (Reference  4)  and  later  modified  by 
Boyer  and  Russell  (Reference  5).  The  program  is  written  in  FORTRAN  and 
consists  of  a  driver  program  and  several  subroutines.  Each  subroutine 
calculates  one  of  the  three  temperature  types  (surface,  interior,  or 
interface  node) ,  or  performs  other  tasks  required  by  the  program.  The  model 
has  been  executed  on  a  16-blt  desk-top  computer.  The  program  is 
interactive;  it  requests  the  user  to  enter  the  data  and  output  file  names. 

Tne  input  data  consist  of  dimensions,  thermal  properties  of  each  of  the 
layers,  flame,  and  cooling  parameters,  and  program  control  parameters.  The 
input  data  required  by  the  program  are  listed  in  Appendix  C. 

B.  COMPARISON  WITH  OTHER  MODELING  AND  EXPERIMENTAL  DATA 

The  model  was  tested  by  executing  several  cases  to  ensure  that  cookoff 
times  were  predicted  correctly.  All  of  the  cookoff  tests  were  conducted 
with  ordnance  engulfed  in  JP-5  fires.  The  model  results  were  first  compared 
to  results  obtained  by  Russell  and  Canfield  (Reference  4)  with  their  flat- 
plate  ordnance  model.  Reference  4,  Appendix  B,  Table  3,  Case  1  was 
used  as  a  baseline.  Russell  and  Canfield  predicted  a  cookoff  time  of  229.7 
seconds.  The  current  model  predicts  a  cookoff  time  of  240.2  seconds. 

Russell  and  Canfield  also  presented  data  for  Mark  82  bombs  with  what  was 
then  the  standard  hot-melt  thickness.  The  average  time  for  reaction  for  the  three 
tests  (Reference  4,  Appendix  B,  Table  5)  was  196  seconds.  Measured  hot-melt  layer 
thicknesses  were  not  available  for  these  tests,  but  radiographs  were  taken 
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of  10  bombs  randomly  selected  from  different  production  runs.  The  average 
of  the  hot-melt  thicknesses  for  the  10  bombs  was  0.33  cm.  Using  the  value 
of  0.33  cm  for  the  hot-melt  thickness,  the  computed  reaction  time  for  a  Mark 
82  bomb  using  the  current  model  is  184.5  seconds. 

C.  W.  Morris  (Reference  8)  also  presented  cookoff  data  that  were  used 
for  comparative  purposes.  Table  1  presents  a  summary  of  experimental  and 
comparative  predicted  cookoff  times  for  five  different  cases. 


TABLE  1.  COOKOFF  TIMES  FOR  MARK  82  AND  MARK  84  BOMBS 
WITH  AND  WITHOUT  FM-26  ABLATIVE  COATING 


r  ~ 

Case 

Ordnance 

FM-26 

thickness , 
mils 

Hot-melt 
thickness , 
mils 

Experimental 
cookoff 
time , 
s 

Morris 
model 
time , 
s 

NMERI 
model 
time , 

5 

1 

Mark  82 

none 

125 

180 

169 

i 

182 

2 

Mark  82 

120 

250 

628 

640 

624 

3 

Mark  84 

none 

125 

182 

196 

1 

228 

4 

Mark  84 

none 

300 

309 

264 

320 

5 

Mark  84 

60 

300 

525 

528 

526 

Agreement  between  the  model  and  experimental  results  is  generally  good.  For 
the  remainder  of  this  report,  ordnance  with  the  FM-26*  ablative  coating  will 
be  referred  to  as  coated  ordance,  and  those  without  will  be  referred  to  as 
uncoated . 


*FM-26  is  an  ablative  coating  manufactured  by  AVCO  Corporation, 
Wilmington,  Massachusetts  01887.  This  material  is  used  as  an  Insulator,  for 
thermal  protection  on  most  nonpropelled  Navy  ordnance. 
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c. 


EFFECTS  OF  COOLING 


When  coolant  is  applied  to  an  engulfed  ordnance,  extension  of  the 
cookoff  time  can  be  expected.  The  cookoff  time  extension  depends  on  many 
parameters.  The  parameters  contributing  most  effectively  to  the  cookoff- 
time  extension  are  the  overall  heat  transfer  coefficient,  cooling  initiation 
time,  and  coolant  temperature. 

A  series  of  calculations  was  made  to  determine  the  dependency  of  the 
cookoff  time  on  the  overall  heat  transfer  coefficient  and  the  start  of 
cooling  time.  The  values  chosen  for  the  overall  cooling  coefficient, 
were  taken  from  the  result  of  the  subscale  calorimeter  testing  reported  in 
sections  IV  and  V.  The  dependency  of  the  cookoff  time  on  the  overall  heat 
transfer  coefficient  and  the  start  of  the  cooling  time  is  shown  in  Figures 
3,  4,  and  5.  Each  of  these  figures  shows  a  composite  cookoff  time 
prediction,  based  on  the  calculation,  for  an  uncoated  ordnance. 

Figure  3  is  a  composite  illustration  of  the  cookoff  time  prediction  for 

an  uncoated  Mark  82  bomb.  The  figure  shows  that,  for  any  U  greater  than 
-5  2 

10  cal/cm  s  K,  ordnance  cookoff  can  be  prevented  if  cooling  starts  at  any 

time  before  about  150  seconds  after  fire  Initiation.  For  the  same  coefficient, 

if  cooling  starts  160  seconds  after  fire  Initiation,  cookoff  can  be  delayed 

approximately  40  seconds.  If  the  cooling  initiation  is  delayed  another  10 

seconds,  i.e.,  until  170  seconds,  the  increase  in  cookoff  time  is  negligible 

(less  than  10  seconds).  A  similar  trend  is  shown  for  the  heat  transfer 

-4  2 

coefficient  of  10  cal/cm  s  K.  Cookoff  can  be  prevented  if  cooling 

starts  at  any  time  before  155  seconds  after  fire  initiation. 

Similarly,  Figures  4  and  5  illustrate  the  cookoff  time  versus  cooling 
initiation  predictions  for  uncoated  Mark  84  general-purpose  bombs  with  different 
hot-melt  thicknesses.  Analogous  to  the  analysis  performed  on  Figure  3,  one  can 
estimate  the  cookoff  time  for  various  values  of  U  .  The  key  to  safe  and 

03 

successful  cookoff  prevention  is  to  apply  the  coolant  before  the  exothermic 
runaway  decomposition  of  the  explosive  inside  the  ordnance  begins. 
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Cookoff  Time  Versus  Cooling  Initiation  Time  Profile 
(Mark  82,  Uncoated,  Table  1,  Case  1). 


350 


Figure  5.  Cookoff  Time  Versus  Cooling  Initiation  Time  Profile 
(Mark  84,  Uncoated,  Table  1,  Case  4). 


SECTION  IV 


SUBSCALE  TESTING 


A .  INTRODUCTION 

Specially  designed  calorimeter  tests  were  conducted  to  obtain  overall 
heat  transfer  coefficients  essential  to  the  cooling  module  of  the  computer 
model.  The  tests  were  also  designed  to  validate  the  performance  of  the 
model.  The  calorimeters  were  placed  in  turbulent  fires  from  a  stagnant  pool 
and  cooled  with  various  cooling  agents. 

Water,  AFFF,  Halon  2402,  and  liquid  nitrogen  were  investigated  as 
cooling  agents.  Water  was  chosen  as  a  baseline  agent  because  of  its  cooling 
capacity  and  its  use  in  previous  cooling  studies.  AFFF  was  chosen  because 
it  provides  a  blanket  which  covers  and  secures  the  fuel  Instead  of  carrying 
it  to  other  areas  as  water  does  (Reference  9).  AFFF  is  also  the  most  likely 
agent  to  be  used  in  a  real  fire  situation.  Studies  show  that  AFFF  does  not 
cool  ns  effectively  as  water.  However,  It  extinguishes  the  fire  while 
cooling  the  ordnance,  and  may  be  a  better  overall  control  agent  for  that 
reason  (Reference  10) .  Halon  2402  was  tested  because  of  Its  superior 
ability  to  provide  rapid  knockdown,  three-dimensional  effectiveness,  and 
slight  cooling  capacity.  It  was  thought  that  the  Halon  2402  would 
extinguish  the  fire  rapidly  while  providing  the  initial  cooling.  Liquid 
nitrogen  was  chosen  as  a  cooling  agent  for  coated  ordnance  because  the 
traditional  cooling  agents  (water  and  AFFF)  were  unable  to  provide  any 
significant  cooling  to  coated  ordnance  during  earlier  test  programs. 

B.  PREVIOUS  WORK 

Experimental  and  theoretical  studies  of  the  response  (cookoff  time)  of 
uncoated  ordnance  when  exposed  to  pool  fires  have  been  reported  by  many 
investigators  (References  9,  10,  11,  12,  13,  and  14).  Specifically,  Hontgas 
(Reference  9)  has  reported  a  series  of  tests  In  which  inert  ordnance  were 
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engulfed  in  pool  fires  and  cooling  was  provided  with  water  and  AFFF  using 
handlines  and  deck  nozzles.  Included  in  the  documentation  are  the 
temperature  history  curves  of  various  locations  within  the  inert  ordnance. 
Multiple  parameters  were  varied  in  each  test  making  it  Impossible  to  get 
correlation  of  the  data  at  different  coolant  flow  rates  or  handline 
locations  relative  to  the  ordnance.  This  study  indicates  that  application 
of  AFFF  would  not  cool  the  ordnance  until  the  fire  had  been  extinguished. 

Cragin,  Pakulak  and  Vernon  (Reference  10)  calculated  a  heat  transfer 
parameter  from  heated  ordnance  to  water  and  AFFF  using  inside  skin 
temperature  of  the  ordnance.  The  procedure  used  during  this  study  was  to 
heat  Inert  ordnance  with  a  propane  burner  rack  to  a  specified  temperature, 
turn  the  burner  rack  off,  and  begin  cooling.  An  extensive  series  of 
photographs  was  Included  in  this  report  to  document  the  effectiveness  of 
the  coolant  coverage  over  the  ordnance.  They  show  that  the  most  effective 
coverage  of  the  ordnance  is  provided  when  the  coolant  is  dispersed  in  a  fog 
pattern.  The  study  shows  that  AFFF  cannot  completely  cover  the  surface  of 
uncoated  ordnance  (especially  the  bottom  of  the  back  side) ,  although  it  can 
cover  coated  ordnance.  This  study  also  concluded  that  if  the  fire  engulfing 
a  coated  Mark  82  is  extinguished  within  7  minutes  the  ordnance  will  not  cook 
off,  even  if  uncooled. 

C.  TEST  SETUP,  INSTRUMENTATION,  AND  MATRIX 

Testing  was  performed  by  suspending  a  calorimeter  0.91  meters  above  a 

2  2 

pool  of  burning  JP-4.  The  pool  size  area  ranged  from  91. A  m  to  274  m  . 

The  calorimeter  consisted  of  three  cylinders  divided  into  five 
circumferential  sections  each.  The  cylinders  and  sections  were  thermally 
insulated  from  each  other  to  allow  quantification  of  the  heat  transfer  rates 
to  various  areas  of  t'..e  calorimeter.  Each  of  the  circumferential  sections 
in  two  of  the  three  cylinders  (the  middle  and  one  end  cylinder)  were 
Instrumented  v;ith  thermocouples  to  measure  the  transient  temperatures 
experienced  during  the  fire  and  subsequent  cooling.  The  third  section  was 
not  Instrumented  because  It  was  symmetric  with  the  first  end  section.  The 
area  with  the  highest  heat  flux  is  the  bottom  of  the  ordnance.  Since  it  Is 
not  always  possible  to  cover  the  bottom  of  the  back  side  (the  near  side  is 
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where  the  coolant  stream  was  directed)  of  the  ordnance  with  coolant,  this 
area  is  of  critical  importance.  This  quadrant  of  the  calorimeter  was  split 
to  have  two  temperature  measurements  In  this  area.  The  construction  of  the 
calorimeters  is  shown  in  Figures  6  through  8. 

Testing  was  first  conducted  using  water,  AFFF,  and  Halon  2402  to  cool 
an  uncoated  calorimeter.  The  calorimeter  was  then  coated  with  0.406  cm  of 
FM-26  and  retested  with  water,  AFFF,  and  liquid  nitrogen.  The  temperatures 
were  read  by  a  datalogger  at  5-second  intervals  and  stored  on  floppy  disks. 

The  data  were  later  reduced  using  a  desk-top  computer.  All  of  the  tests  were 
filmed  with  a  video  camera.  Table  2  shows  the  test  matrix  used  to  evaluate  the 
various  coolants. 


TAELE  2.  COOLING  TEST  MATRIX 


Calorimeter 

Coolant 

Cooling  Flow 

— 

Application 

Type 

Rate 

Method 

uncoated 

water 

341 

L/min 

stream 

uncoated 

^  a 
water 

341 

L/min 

fog 

uncoated 

^  a 
water 

170 

L/min 

fog 

uncoated 

AFFF® 

341 

L/min(water) 

fog 

uncoated 

Halon  2402^^ 

2.5 

kg/s 

stream 

coated 

a 

water 

341 

L/min 

fog 

coated 

^  a 
water 

170 

L/mln 

fog 

coated 

/.FFF® 

341 

L/min(water) 

fog 

coated 

liquid  nitrogen 

g 

The  water  and  AFFF  were  dispensed  using  a  variable-cone,  variable-flow¬ 
rate  nozzle  attached  to  a  3.81-cir.  handline. 

'^The  Halon  2402  was  dispensed  from  a  pressurized  system  through  a  standard 
USAF  0.71-cm  halon  nozzle. 

Q 

During  the  liquid  nitrogen  tests  the  Dewar  was  pressurized  to  1000  kPa. 
The  liquid  discharge  valve  was  then  opened  and  the  nitrogen  flowed  at  whatever 
rate  the  Dewar  could  maintain.  By  the  time  the  nitrogen  reached  the  diocharge 
of  the  feed  pipe  (1.83  meters  long  by  2.54  cm  diameter),  almost  all  of  it  had 
vaporized  and  the  calorimeter  was  engulfed  In  vapor. 
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Figure  6.  Calorimeter  Engineering  Drawing. 


(a)  Large-Scale  Calorimeter  on  Test-Bed. 


b)  Fire,  Coolant,  and  Calorimeter  Interaction 


Figure  8.  Calorimeter  Assembly  and  Testing. 


D.  OVERALL  HEAT  TRANSFER  COEFFICIENT 


The  overall  heat  transfer  coefficient  can  be  calculated  from  equating 
the  convective  heat  transfer  with  the  conductive  heat  losses  at  the 
calorimeter  surface; 


U  (T  -  T  )  -  K  dT 
oa  w  s  ~ 

dr 


In  Equation  (16)  the  variable  U  represents  the  overall  heat  transfer 

oa 

coefficient,  and  T  and  T  arc  coolant  and  surface  temperature, 

^  ^  dT 

respectively.  In  the  right-hand  side,  K  is  the  thermal  conductivity  and  ~ 

is  the  spatial  temperature  gradient  across  the  calorimeter  wall  and  can  be 

calculated  from  the  Fourier  Conduction  Law.  However,  the  total  heat  loss 

from  the  calorimeter  during  the  cooling  period  is  given  by 


q  -  ^  (17) 

^  at 

where  p  and  C  are  density  and  the  heat  capacity  respectively,  which  are 

^  aT 

known  quantities  for  the  calorimeter.  The  term  ^  is  the  slope  of  the 

a  t 

temperature  history  profile,  which  is  known  at  any  time.  Since  all  parts  of 
Equation  (17)  are  known,  the  heat  loss  -ate,  q,  can  be  calculated. 
Substituting  q  from  Equation  (17)  into  the  Fourier  Conduction  Law,  one 
obtains 


q/A  -  K  ^ 
dr 


The  left-hand  side  of  Equation  (18)  is  known  and  Equation  (16)  can  be 
writter  as 


U  (T  -  T  )  -  K  dT  -  q/A 
oa  w  s  T* 

dr 
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A  computer  program  was  written  to  calculate  from  the  temperature 
history  profiles  and  is  presented  in  Appendix  C. 

Figure  9  illustrates  the  calorimeter's  cross  section  and  the  location 
of  thermocouples.  The  arrows  specify  the  direction  in  which  the  coolant  was 
applied.  F.ffects  of  various  coolants  were  tested  experimentally.  Figures 
10  and  11  show  typical  temperature/tlme  profiles  for  heating  and  cooling. 

In  this  case  the  cooling  agent  is  water.  The  temperature  profiles  resulting 
from  other  coolants  are  given  in  Appendix  D.  Discussion  of  each  coolant 
follows. 

The  test  results  showed  that  of  the  two  water  flow  rates  tested,  the 

higher  the  water  flow  rate  the  more  effective  the  cooling.  The  overall 

heat  transfer  coefficients  for  the  different  segments  of  the  calorimeter 

-2  -3 

for  a  water  flow  rate  of  341  L/min  ranged  from  10  to  10  cal/cm*'  s  K, 

-2  -5  2 

while  they  ranged  from  10  to  10  cal/cm  s  K  for  a  water  flow  rate  of 
170  L/min.  This  is  reflected  in  the  temperature  profiles.  The  temperature 
profiles  for  the  341  L/min  flow  rate  show  an  immediate  and  dramatic 
temperature  reduction  for  all  parts  of  the  calorimeter  (Figures  D-1  and  D-2). 
The  temperature  profiles  for  the  170  L/min  flow  rate  show  rapid  and  dramatic 
temperature  reductions  only  at  the  point  where  the  water  struck  the 
calorimeter  (Figures  10,  11,  D-3,  and  D-4) .  There  was  a  delay  in  the  onset 
of  the  temperature  drop  in  other  parts  of  the  calorimeter  and  the  temperature 
drop  rate  was  not  nearly  as  high  as  with  the  341  L/min  case.  The  cooling 
was  ineffective  on  some  parts  of  the  calorimeter. 

AFFF  rapidly  extinguished  the  fires  during  testing.  It  then  cooled 
the  calorimeters,  although  at  a  rate  much  lower  than  that  of  water  (Figures 
D-5  through  D-8) .  The  AFFF  was  incapable  of  providing  cooling  to  some  parts 
of  the  calorimeter.  The  overall  heat  transfer  coefficients  for  AFFF  at 
341  L/min  ranged  from  10  to  10  cal/cm  s  K.  The  overall  heat  transfer 
coefficients  for  170  L/min  were  similar. 
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Figure  11.  Calorimeter  Temperature/Time  Profile  (Uncoated,  Water  Coolant,  170  L/min) . 


The  Halon  2402  was  ineffective  as  a  coolant.  It  was  able  to  stop  the 

temperature  rise  at  the  points  where  it  struck  the  calorimeter,  but  was 

unable  to  cool  other  areas  (Figures  D-9  and  D-10) .  The  overall  heat  transfer 

-4  2 

coefficients  ranged  from  0  to  10  cal/cm  s  K.  The  Halon  2402  did  not 
extinguish  the  fire  since  was  net  sprayed  at  the  base. 

The  testing  on  coated  calorimeters  showed  the  remarkable  insulating 
ability  of  the  FM-26  coating.  Before  cooling  was  initiated,  the  temperature 
rise  was  a  fraction  of  that  for  the  uncoated  calorimeter  (Figures  D-11 
through  D-14).  The  insulating  ability  of  the  FM-26  also  dramatically  slows 
down  the  rate  of  energy  removal.  The  overall  heat  transfer  coefficient  with 
a  water  coolant  flow  rate  of  341  L/min  ranged  from  10  to  10  cal/cm  s  K. 

The  liquid  nitrogen  provided  very  effective  cooling  in  the  area  where 

the  nitrogen  impinged  on  the  calorimeter  (Figure  D-15).  The  fixed-pipe 

feed  system  for  the  liquid  nitrogen  did  not  allow  it  to  be  spread  over  the 

calorimeter  surface,  although  the  nitrogen  did  provide  more  effective  cooling 

than  the  water  on  the  impinged  area.  The  overall  heat  transfer  coefficients 

-2  -5  2 

for  the  liquid  nitrogen  ranged  from  10  to  10  cal/cm  s  K. 
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SECTION  V 


CONCLUSIONS  AND  RECOMMENDATIONS 

Response  of  ordnance  to  accidental  flame  Impingemant  was  studied 
experimentally  as  well  as  numerically.  Calculations  were  made  to  quantify 
the  effect  of  the  various  coolants  on  cookoff  times.  Experiments  were 
conducted  to  assess  the  fire  interaction  with  ordnance  to  quantify  effective 
coolant  application  rates,  to  evaluate  various  cooling  agents,  and  to 
provide  baseline  heat  transfer  coefficients  essential  to  cookoff  prediction 
calculations . 

Calculations  show  that  extension  of  the  cookoff  time  can  be  obtained 
when  coolant  is  applied  to  the  ordnance  as  long  as  the  temperature  of  the 
explosive  has  not  reached  the  critical  point,  that  Is,  the  point  where  an 
irreversible  exothermic  reaction  begins.  It  is  concluded  that  cookoff 
time  extension  strongly  depends  on  the  overall  heat  transfer  coefficient 
and  the  cooling  initation  time. 

The  effects  of  various  coolants  on  cookoff  time  were  evaluated 
experimentally.  Water,  AFFF,  Halon  2402,  and  liquid  nitrogen  were 
investigated  as  cooling  agents.  Although  AFFF  does  not  cool  as  effectively 
as  water,  it  extinguishes  the  fire  while  cooling  the  ordnance  and  may  be 
considered  a  superior  agent.  Ordnance  with  intumescent  coatings  were 
compared  with  uncoated  ordnance.  Limited  experimental  data  suggest  that 
FM-26  coating  is  effective  initially  in  preventing  heat  from  transferring 
in,  but,  at  the  same  time,  the  char  layer  formed  when  the  coating  is  heated 
hinders  removal  of  heat  by  the  coolant  and  does  not  significantly  alter  the 
cookoff  time. 

It  is  recommended  that  additional  studies  be  conducted  to  further 
evaluate  the  effectiveness  of  water  and  AFFF  on  cookoff  time  of  coated 
ordnance  in  terms  of  the  physical  properties  of  various  coating  materials. 
This  would  require  a  combination  of  predictions  using  i.ie  computer  codes 
developed  and  carefully  instrumented  tests  using  live  ordnance. 
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program  mad 
$scorage : 2 
$nofloatcalls 

real  hm ,hml ,hm2 

character* 10  matl ( 5) , explid, weapid ,title(8) 
common  /  a  /  t(1000. 2)  ,  tm(5)  ,  qtn(lOOO)  ,hm(5)  ,hrrl(5)  ,hm2(5)  ,k, 

&  kk(5) .ll,12,13,m.l,pt(5) .ipt(5).rad(1000) ,d(5) .keyeqn. Icrit . 

&  c tr , maxt . dx ( 5) , dam( S) , n , Ipt , weapid , explid, matl , inc , r , title , 

&  V, time . tO, j J J (16) 

common  /  tp  /  ckl , ck2 , ck3 . fOl , f02 , f03 , rfk(5) , rcx( 5) , alpha(5) , 

6.  c(5) .cone(5) .cchree(5) , first(5) ,second(5) , thlrd(5) ,fourth(5) , 

&  sevnth(5).eighth(5).ctwo(5) .fifth(5),8lxth(5) ,tmld(5) .thlgh(5) , 

&  ae,rho(5)  ,xr.(5,5)  ,q(5,5)  ,e(5,5),z(5.5)  .keytherra,  theta(5)  .ck(5) 
common  /  be  /  tf lame , t£ , epsilon. sigma, tc , tlm, dt, tbld, he , rt ,bi . 

6i  thlghst(20) , layrout, Irestart , tcooi , timcool ,uoa, tw 
tread  -  0. 

10  call  readtnpcC tread) 
tf  (  k  .eq.  0  )  go  to  80 
tf  (  Irestart  .ne.  1  )  go  to  20 
go  to  30 
20  call  fourier 
call  Inltllze 

c . evaluate  temperatures  In  layers  and  interfaces . 

Ipt  -  1 

30  do  60  mj-lpt,maxt 
1-2 

tlm  -tim+dt 
m-0 

call  layer 

tf  (Icrit  .eq.  1)  call  prntcrlt 
if  (icrit  .eq.  1)  go  to  70 

c . check  for  printout  time . - 

if (amod(ctr ,v) )  50,40,50 
40  call  prntrslt 
50  ctr-ctr+1. 

do  60  iset-l,n 
t ( Ise t ,l)-t(iset,2) 

60  continue 

70  go  Co  10 

80  stop 

end 

subroutine  readlnpt( iread) 
characcer*64  file?,  flle6,  file9 

charac ter*10  matl ( 5) , explid, weapid ,tltle(8) 
common  /  a  /  t(1000,2) , tm(5) ,qm(1000) ,hm(5) ,hml(5) ,hm2(5) ,k, 

&  kk(5) , 11, 12 ,13 ,m, l.pt(5) ,ipt(5) ,rad(1000) ,d(5) .keyeqn, icrit, 

&  ccr , maxt , dx(5) , dam(5) , n, Ipt .weapid , explid, matl , inc , r , title , 

&  V, time , tO , j j J (16 ) 

common  /  tp  /  ckl , ck2 , ck3 , fOl , f02 , f03 , rfk(5) , rcx( 5) , alpha( 5) , 

6c  c(5) . cone (5) ,cthree(5) , f Irs t ( 5) , second ( 5) , third(5) . fourth(5) . 

&  sevncb(5) , eighth (5) ,ccwo(5) , fifth(5) , sixth (5) , tmid(5) , chigh(5) , 
6.  ae ,rho(5) ,xm(5,5),q(5,5),e(5,5),z(5,5) , key therm, theta (5) ,ck(5) 
common  /  be  /  tf lame , tf , epsilon, sigma, tc , tim , dt , tbld, he , rt ,bi , 

&  thighst(20), lay  rout, ire  start,  icool ,  tirncool .  uoa,  tw 
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c . read  in  computation  requirements . 

c  k  -  number  of  layers,  v  -  calc . /pj^lntout 

c  time  —  max.  time  for  calc, 

c  dt  -  time  increment 


if  (  iread  .ne,  0.  )  go  to  200 

iread  -  1. 

write(*,100) 

100  formate*  Enter  name  of  input  data  file:  '\) 
read(*. ' (a) ' )  file? 
write(’f.llO) 

110  formate*  Enter  name  of  file  for  output  data:  '\) 
reade*.‘ea)')  flle6 
writee*, 120) 

120  formate*  Enter  name  of  restart  data  file:  '\) 
reade*, ‘ ea) * )  file9 
opane?, file-file? , status-*old* ) 
openefi 1 f ile-f ile6 , status— *  new* ) 

200  reade?.*)  k.  r 

if  ek.eq.O)  return 
read  e^.*)  irestart , icool , layrout 
if  e  Irestart  .ne.  1  )  go  to  5 
opener, file-file9, status—' old* ) 
go  to  8 

5  opene9,flle-flle9 .status- 'new* ) 

8  reade?, 1)  etitle( j ) . j-1.8) 

reade?,*)  v,  time,  dt,  inc 
do  10  j  -  l.k 
reade?, 1)  mitiej) 

10  reade?,*)  dame j ) , rhoe j ) . cej ) ,kkej ) , pte j ) ,hme J ) , tme J ) 

do  20  J  -  1.  k 

raade? . *)  coneej ) , f irstej ) , seconde j ) , third ej ) , fourthej ) 
reade?,*)  ctwoej),  fifthej),  sixthej ) , tmidej ) 

reade? ,*)  cthreee j ) . sevnthej ) .elghthej ) , thighej) 

20  continue 

reade?,*)  tO.tc.epsilon.di, t flame , tbld 
rt  -  di  /  2. 
reade?, 1)  explid 
reade?,  1)  weapid 
do  30  ki  -  l.k 

J  -  1 

reade?,*)  eeki , j ) , z eki , j ) , xmeki , J ) , qekt , j ) 
if  ekk(kl)  .eq.  0  )  go  to  30 
do  25  j  -  2,kkekt) 

25  reade?,*)  e  eki  ,j  ),  z  eki  ,j  ),  xmeki  .  D.  qeki  ,J  ) 

30  continue 

if  e  icool  .ne.  1)  go  to  UO 
reade?,*)  timcool , uoa, tw 
^0  lf(  irestart  .ne.  1  )  return 
reade9,2)  n, tim.maxt.lpt 
rcade9.3)  etei  ,1) , 1-l.n) 

1  format  eSalO) 

2  formateih  , ilO , 5x , f ? . 3 , 5x , ilO, 5x , ilO) 

3  formateih  ,8eix,f9.5)) 

return 


A3 


end 

subroutine  initlize 

character*10  matl(5) .explid.weapid, tltle(8) 
conunon  /  a  /  t(1000 , 2) ,  tm(5)  ,qm(1000)  ,hm(5)  .hml(5)  .hni2(5)  ,k, 

&  kk(5)  ,11 , 12,13  .tn.  i  ,pr.(5)  ,ipt(5)  ,rad(1000)  ,d(5)  .keyeqn,  icrit, 

&  ctr ,maxc ,dx(5) ,dam(5) . n , Ipt , weapid , explld ,matl , tnc . r , 1 1  tie , 

&  v,uime,t0,jjj(16) 

common  /  tp  /  ckl , ck2 , ck3 , fOl , f02 . f03 . r fk( 5) , rcx( 5) , alpha( 5 ) , 

&  c(3)  cone(5),cthree(5),first(5),second(5),third(5).fourth(5) , 

&  sevntb(5) ,eighth(5) .ctwo (5) , f if th(5) . slxLh(5) . tmld(5) , thigh( 3) . 

&  ae,rho(5) ,xm(5,5) ,q(5,5) . e(5 . 5) , z (5 . 5) , keytherm, theta(5) ,ck( 5) 

common  /  be  /  tf lame , tf , epsilon, sigma, tc , tim.dt, tbld ,hc , rt ,bl , 

L  thighEt(20) , layrout ,  irestart . tcool , timcool ,uoa, tw 
icrit  -0 
etr  -0.0 

sigma  -  1.356e-12 
maxt  -if ix(time/dt) 
tim-- dt 

13  -  k-1 

do  10  ml-1.13  , 

10  rcx(ml)-2 .*dt/(rho(ml)'*c(ml)*dx(ml)+rho(ml+l)*c(ml+l)*dx(ml+l) ) 

do  20  kl— l,k 

hml(kl)-rho(kl)*dx(kl)*. 5*hm(kl) 

20  hm2(kl)  -  2.*hml(kl) 

c . evaluate  radius  at  every  point . 

rad(l)  -  rt 
13-1 . 

16-0. 

do  160  iA-l,k 
i6-ipt(i4)+16 
do  165  mp-i5 , 16 

165  rad(mp+l)-rad(mp) -dx(i4) 

160  15-16+1 

c . set  temperatures  initially  to  a  constant . 

do  201  init-l,n 
qm( init)— 0 . 
t ( ini t , 2 )-c0 
201  t(init,l)-t0 

do  1609  ii-1,20 
1609  chighst(ii)  -  0.0 
return 
end 

subroutine  fourier 

character *10  matl(5), expl id , weapid , t i tl e  (  8 ) 
common  /  a  /  t  ( 1000 , 2 )  ,  tm(  5)  ,  qmUOOO)  , hm ( 5)  , hml  (  5  )  ,  hm2  ( 5)  ,  k  , 

6.  kk(5)  ,ll,12,13,m,i,pt(5),ipt(5), r ad ( 1000) , d( 5) , keyeqn , Ici it , 

&  ctr , maxt .dx(5) ,dam(5) , n, Ipt , weap id , expl id , mar 1 , inc , r , title , 

6  v,time,t0,jjj(16) 

common  /  tp  /  ckl , ck2 , ckl , f 01 , f 02 , f03 , r f k( 5 ) , rex { 5 ) , alpha ( 5 ) , 

&  c (5) , cone (b),cthree(5),firsc(5), second( 5) , third(5) , fourth ( 5) , 

&  sevnLh(5) ,eighth(5) , c two ( 5) . f if th ( 5 ) , s ixth( 5 ) , tmid( 5 ) ,thigh(5) , 
&  ae,rho(5) ,xm(5, 5) ,q(5,5),e(5,5).z(5,5) , key therm , the ta ( 5) ,ck(5) 
common  /  be  /  t flame , tf , eps lion , sigma , tc , tim , dt , tbld , he , rt , bi , 

&  thighst (20) , layrout .irestart, icool , timcool ,uoa, tw 
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c . evaluate  and  stablize  fourier  modulus . 

45  tl  -  t0*1.8  -460. 

he  -  0.00134*(21.666/(rt*2.))**.195 
n-l .  0 

do  15  j-l,k 

ck(J )-(cone(j )+f irst(J )*tl+3econd( j )*tl*tl+Chird( J )*tl*tl*tl+ 
IfourthCJ )*(tl**4) )/241.9 
alpha(J )-ck(j )/(rho(J )*c(j ) ) 
tpt(j)  -  ifix(pt(J)) 

42  dx(J)-dam(j)/pt(j) 
tf(dt.ne.0)go  to  27 
do  28  ij-l.k 
lf(lj .eq.l)go  to  41 

ck( ij )- (cone(lJ )+f lrst(tj )*tl+8econd( IJ )*tlwtl+thlrd( IJ )*tl*tl*tl 
l+fourth(lJ)*(tl**4))/241.9 
alpha ( 1 J ) -ck ( I J ) / ( rho ( 1 j ) *c ( i J ) ) 
ipt(tj)-lflx'pt(ij)) 
dx(ij)-dam  i  '/pt(lj) 
d(lj )-0 . 2*(dx(lJ )**2)/alpha(lJ ) 
go  to  39 

41  bl-hc*dx(l)/ck(l) 

d(tj)-0.5*(dx(l)**2)/((l.+bl)*alpha(l)) 

39  lf(d(l) .le.d(lj))go  to  28 
d(l)-d(lj) 

28  continue 
dt-d(l) 

27  theta(J )-alpha(j ) *dt/(dx(J )**2) 
rfk(J )-theta(j)/ck(J) 

If ( j . eq . 1)  go  to  20 
go  to  21 

20  bl-hc*dx(l)/ck(l) 
if(theta(l)-0.5/(l.+bl))  112,112,22 

21  if ( theta( J ) . le . 0 . 2 )  go  to  112 

22  theta( j )-0 . 75*theta( j ) 

dx(J )-sqrt(alpha(J )*dt/theta(J ) ) 
ipt(J)-lflx(dam(J)/dx(J)) 
pt(j)-float(ipt(J)) 
go  to  42 
112  n-n+ipt(J) 

15  continue 

return 

end 

subroutine  layer 

character^lO  macl(5) , explld,weapid , tltle(6) 
common  /  a  /  t(1000, 2) , tm(5) ,qm<1000) ,hm(5) ,hml(5) ,hm2(5) ,k, 

&  kk(5) ,ll,12,13,m,l,pt(5),ipt(5), rad (1000) , d(5) ,keyeqn, Icrlt , 

&  ctr ,maxt , dx(5) , daro(S) ,n, ipt ,weapld,explld,matl , Inc , r , title , 

&  v,tlme,t0,jjj(16) 

common  /  tp  /  ckl , ck2 , ck3 , fOl , f02 , f03 , rfk(5) , rcx(5) , alpha(5) , 

6.  c(5) , cone (5) ,cthree(5) , first (5) ,second(5) ,thlrd(5) ,fourth(5) , 

&  sevnth(5) ,elghth(5) ,ctwo(5) ,flfth(5) ,8lxth(5) ,tmld(5) ,thlgh(5) , 

&  ae  ,  rh<j(5)  ,xm(5 , 5)  ,q(5,5),e(5,5),z(5,5)  ,keytherra,  theta(5)  ,ck(5) 
common  /  be  /  tf lame , tf , epsilon, s Igma , tc , tlm, dt , tbld , he , rt ,bl , 

&  thlghst(20), lay rout , Ires tart , icool , tlmcool ,uoa , tw 
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do  100  11-1,  k 
12-11+1 

if  (  11  .gt.  1  )  go  to  10 
call  tsurf 
10  m  -  m  +  ipt(ll) 
if  (  11  .eq,  1  )  go  to  20 
i  -  1  +  2 
20  call  thermprp 
keyeqn  -  0 

if  (  t(l.l)  .It.  tm(ll)  )  keyeqn  -  1 
if  (  qm(i)  .ge.  hm2(ll)  )  keyeqn  -  1 


call  newtemp 
dtmin  -  f02  *  0.01 

if  (  abs(  t(i,2)  -  t(i.l)  )  .le.  dtmin  )  go  to  60 
30  if  (  1  .eq.  m)  go  to  40 
i  -  i  +  1 


go  to  20 

40  if(  11  .eq.  k  )  go  to  50 
call  tintface 


CO  to  60 

c  ***  calculate  temperature  of  center  of  bomb  +********************** 
50  1  -  1  +  1 

t(i,2)  -  f02  *  (  (  dx(ll)**2  )  *  ae  /  ck2  +  2.  +  t(l-l,l)  + 

’&(!./  f02  -  2.  )  *  t(i.l)  ) 

c  ***  check  to  see  if  bomb  has  cooked  off  ************+************** 
60  if  (  11  .ne.  k)  go  to  100 
m2  -  m  •  ipt(ll)  +  1 
do  70  11  -  m2,m 

if  (  t(il,2)  .It.  tc  )  go  to  70 
icrlt  -  1 
return 
70  continue 
100  continue 


return 

end 

subroutine  tsurf 

characcer*10  matl(5) , explid .weapld, title(8) 
common  /  a  /  t(1000, 2) . tm(5) ,qm(l000) ,hm(5) ,hml(5) ,hm2(5) ,k, 

(.  kk(5)  ,11, 12,13,  m,l,pt(  5),  lpt(5).ra<l(1000),d(5), keyeqn,  icrlt, 

&  ctr  ,maxt ,  dx(5)  ,  dani(5)  ,n,lpt , weapld,  explid, mat  1 ,  Inc ,  r,  title  , 

&  v,time,t0, jjj(16) 

common  /  tp  /  ckl , ck2 , ck3 . fOl , f02 . f03 . rfk(5) , rcx( 5) , alpha(5) . 

&  c(5) , cone (5) ,cthree(5)  .first (5) .second(5) ,third(5) ,fourth(5) , 

&  sevnth( 5) , e Ighth ( 5) , c two( 5) .fifth (5) .sixth (5) ,tmld(5) .thigh (5), 

&  ae,rho(5) ,xm(5,5) ,q( 5 , 5) ,e ( 5 , 5) , z(5 , 5) .keytherm , theta( 5) , ck( 5) 

common  /  be  /  t flame , tf , eps lion , s igma , tc , tlm , dt , tbld , he , rt , b 1 , 

&  thighst(20) , layrout , Irestart , icool , tlmcool ,uoa , tw 
f(x,a,b.c,d.e)  -  (  a  +  b*-  +  c*x*x  +  d*x*x*x  +  e-*(x**4)  )  /241.9 

g(x',a,b,c,d)  -  <  a  +  b+(  1.  -  oxp(<c-x)/d)  )  )  /  241.9 

zz(x,a,b,c)  -  (  a  +  b/(x*x)  +  c/(x*+4)  )/  241.9 
c  ***  calculate  Internal  energy  generated  at  surface  ****-*************-^* 
if  (  tim  . ge .  tbld  )  go  to  5 

tf  -  (  tflame  -  tO  )  *  (  tlm**3  /  tbld**3  )  +  tO 

go  tc  7 
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5  tf  -  tflame 

7  tavg  -  (  t(l,l)  +  t(2,l)  )  /  2. 
ae  -  0 

do  10  ti  -  l,kk(ll) 

10  ae  -  ae  +  xmdl.ii)  *  q(ll,li)  *  2(11,11)  * 

&  exp(-e(ll,lt))/(r*tavg) 
ae  -  ae  *  rho(ll) 

c  ***  calculate  thermal  properties  of  first  two  nodes  ***************** 
tlf  -  t(l,l)  *  1.8  •  ^60. 
t2f  -  t(2.1)  *  1.8  -  460. 

c  ***  account  for  reaction  of  ablative  coating  ************************ 
lf(  layrout  .ne.  1  )  go  to  14 
If  (  tlf  .  It.  thlghst(l)  )  go  to  11 
thighst(l)  -  tlf 
go  to  12 

11  tlf  -  thlghst(l) 

12  if  (t2f  .It.  thighst(2)  )  go  to  13 
thlghst(2)  -  t2f 

go  to  14 

13  t2f  -thlghst(2) 

14  avgts  -  (  tlf  +  t2f  )  /  2. 

If  (  avgts  .ge.  tmld(ll)  )  go  to  16 

cka  -  f (avgts , cone (11) , first (11) , second (11) , third (11) , fourth (11) ) 
go  to  20 

16  if  (  avgts  .ge.  thlgh(l)  )  go  to  18 
cka  -  2z(  avgts,  ctwo(ll),  flfth(ll) ,  sixth(ll)  ) 
go  to  20 

18  cka  -  g(avgt8,cthree(ll) ,8evnth(ll) ,thlgh(ll) ,elghth(ll)) 
c  ***  branch  to  proper  temperature  calculation  ********************** 

20  If  (  t(l.l)  .It.  tm(l)  )  go  to  30 
If  (  qm(l)  .ge.  hml(l)  )  go  to  30 

c  ***  melting  temperature  and  energy  calculation  ******************** 

If  (  Icool  .ne.  1  .or.  tlm  .It.  tlmcool  )  go  to  25 

qm(l)  -  qm(l)  +  dt  *  dx(l)  *  ae  /  2 .  +dt  *  (uoa  *  (tw  •  t(l,l))* 

2  (rt  dx(l)  /  2.)  *  cka  ★  (  c(l.l)  -  t(2,l)  )  /  (  dx(l)  *  rt  )  ) 

go  to  28 

25  qm(l)  -  qm(l)  +  dt  *  dx(l)  *  ae  /  2 .  +dt  *  (  he  *  (tf  -  t(l,l))+ 

1  epsilon  *  sigma  *  (  tf**4  •  t(l,l)**4  )  - 

2  (rt  .  dx(l)  /  2.)  *  cka  *  (  t(l.l)  -  t(2,l)  )  /  (  dx(l)  *  rt  )  ) 

28  tf  (  qm(l)  .gt.  hml(l)  )  go  to  40 

if  (  qm(l)  .It,  0.0  )  go  to  50 

c(l,2)  -  t(l.l) 

return 

c  ***  calculate  new  temp  rise  from  external  heat  fluxes  ************* 

30  if  (  icool  . ne .  1  .or.  tim  .It.  tlmcool  )  go  to  35 
t(l,2)  -  t(l,l)  +  rfk(l)  *  (  dx(l)**2  )  *  ae  + 

&  2,  *  rfk(l)  *  (  rt  *  dx(l)  *  (  uoa  *  (  tw  -  t(l,l)  )  )  • 

6.  (  rt  -  dx(l)  /  2.)  *  cka  *  (  c(l.l)  -  t(2,l)  ))/  (  rt  -  dx(l)/4.) 

return 

35  t(l,2)  -  t(l,l)  +  rfk(l)  *  (  dx(l)**2  )  *  ae  + 

6.  2.  *  rfk(l)  *  (  rt  *  dx(l)  *  (  he  *  (  tf  -  t(l,l)  )  + 

&  epsilon  *  sigma  *  (  tf**4  -  t(l,l)  **  4  )  )  - 

&  (  rt  -  dx(l)  /  2.)  *  cka  *  (  c(l.l)  -  t(2,l)  ))/  (  rt  -  dx(l)/4.) 
return 
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c  ***  calculate  temp  rise  from  excess  energy  above  melting  energy  *** 

40  t(l,2)  -  t(l,l)  +  2.  *  (  qm(l)  -  hn\l(l))  /  (rho(l)  *  c(l)*dx(l)) 
C!ni(l)  -  hm(l) 
return 

c***  calculate  temp  drop  from  endothermic  reaction  ****************** 

50  c(l,2)  -  t(l,l)  +  2.  *  qm(l)  /  (  rho(l)  ♦  c(l)  *  dx(l)  ) 
qm(l)  -0.0 
return 
end 

subroutine  thermprp 

character*10  macl(5) , explid .weapid, title (8) 
common  /  a  /  t(1000 , 2) , tm( 5) ,qm( 1000) ,hm(5) ,hml(5) ,hm2(5)  ,k, 

&  kk(5) ,11.12.13.m,l,pC(5).lpt(5), rad(lOOO) , d( 5) .keyeqn, icrit , 

&  ctr , maxt , dx( 5) , dam(5)  ,n, Ipt .weapid, expl id ,matl , Inc , r , title , 

&  V, time , tO . J J j (16) 

common  /  cp  /  ckl , ck2 , ck3 , fOl, f02 , f03 , rfk(5) , rcx( 5) ,alpha(5) , 

&  c(5) , cone(5) ,cthree(5),first(5), second(5) , third(5) , fourth(5) , 

6  sevnth( 5) , elghch( 5) . ctwo( 5) .fifth (5) .sixth (5) , traid(5) , thigh (5) . 

&  ac , rho(5) .xm(5,5),q(5,5).e(5.5),z(5.5) . keytherm . the ta(5) . ck(5) 
common  /  be  /  tflame , tf , epsilon , sigma . tc . tim. dt , tbld .he , rt ,bi . 

&  thighst ( 20) , lay rout , ires tart . icool . timcool . uoa. tw 
f (x . a .b , c . d . e)  -  (  a  +  b*x  +  c*x*x  +  d*x*x*x  +  e*(x**4)  )  /241.9 
g(x.a.b,c.d)  -  (  a  +  b*(  1.  •  exp((c-x)/d)  )  )  /  241.9 
zz(x,a,b,c)  -  (  a  +  b/(x*x)  +  c/(x**4)  )  /  241.9 
c  ***  this  .subroutine  calculates  thermal  conductivities  (ck), 
c  ★★★  fourler  moduli!  (fox)  for  each  node  point  and 

c  ***  internal  energy  generation  (ae)  for  the  current  c/ilculation  node. 
If  (  keytherm  .ne.  0  )  go  to  50 
tempi  -  t(l*l,l)  *  1.8  460. 

Cemp2  •  t(l,l)  *  1.8  •  460. 
if  (  11  .ne.  1  .or.  layrout  .ne.  1  )  go  to  6 
If  (  tempi  .le.  thighst(l-l)  )  go  to  4 
thighsc(l-l)  -  tempi 
go  to  6 

4  tempi  -  chighsc(i-l) 

6  if  (  tempi  .ge.  tmid(ll)  )  go  to  8 
ckl  -f( tempi  .cone (11)  .firstdl)  .8econd(ll)  .  third (11)  ,fourth(ll)) 
go  to  20 

8  if  (  tempi  . ge .  rhigh(ll)  )  go  to  10 
ckl  -  zz  (  tempi,  ctwo(ll),  fifth(ll) ,  slxth(ll)  ) 
go  to  20 

10  ckl  -  g( tempi , cthree (11) , gevnth(ll) , thigh(ll) . elghth(ll) ) 

20  fOl  -  rfk(ll)  *  ckl 

if  (  11  .ne.  1  .or.  layrout  .ne.  1  )  go  to  23 
if  (  temp2  . le .  thlghst(l)  )  go  to  22 
thlgh8t(i)  -  temp2 
go  to  23 

22  temp2  -  thtghst(i) 

23  if  (  temp2  . ge .  tmid(ll)  )  go  to  25 

ck2  -f  ( temp2  ,  cone  (ll),flr8t(ll),  second(ll)  .  tVilrd(  11) ,  fourth(  11 ) ) 
go  to  40 

25  if  (  temp2  .ge.  thigh(ll)  )  go  to  30 
ck2  -  zz(  temp2,  ctwo(ll)  ,  fifth(ll) .  slxth(ll)  ) 
go  to  40 


30  ck2  -  g( tempZ , cthree(ll) , sevnch(ll) , thigh(ll) , eighth(ll) ) 

40  f02  -  rfk(ll)  *  ck2 
go  to  60 

50  ckl  -  ck2 
ck2  -  ck3 
fOl  -  f02 
f02  -  f03 

60  ae  -  0. 
do  70  11  -  l.kk(ll) 

70  ae  -  ae  +  xm(ll,il)  *  q(ll,tl)  *  z(ll.ii)  * 

&  exp(-e(ll.ll)/(r*t(l.l))) 
ae  -  ae  *  rho(ll) 
temp3  -  t(t+l,l)  *  1.8  -  460. 

If  (  11  .ne.  1  .or.  layrouc  .ne.  1  )  go  to  76 
If  (  temp 3  .le.  thighs t( 1+1)  )  go  to  74 
thlgh8t(l+L)  -  terap3 
go  to  76 

74  temp3  -  thighst(i+l) 

76  If  (  teiiip3  .ge.  tinld(ll)  )  go  to  78 

ck3  -f (tenip3 , cone(ll) ,flrat(ll) ,second(ll) ,thlrd(ll) ,fourth(ll)) 
go  to  90 

78  If  (  ten>p3  .ge.  thlgh(ll)  )  go  to  80 
ck3  -  zz(  temp3,  ctwo<ll),  flfth<ll),  8txth(ll)  ) 
go  to  90 

80  ck3  -  g(tcmp3 ,cthree(ll) ,savnth(ll) ,thlgh(ll) ,elghth(ll)) 

90  £03  -  rfk(ll)  ★  ck3 
return 
end 

subroutine  newtemp 

character*10  matl(5) ,explid,weapld, tltle(8) 
common  /  a  /  t (1000 , 2) , tm(5) .qm(ioOO) ,hm(5) ,hml (5) ,hm2(5) ,k, 

&  kk(5),ll,12,13,m,l,pt(5),lpt(5). rad (1000) ,d(5) .keyeqn, Icrlt , 

4  ctr ,maxt , dx(5) , dam(5) ,n, Ipt ,weapld,explid,rratl , Inc ,r , title , 

&  V, time , tO , J J J ( 16) 

common  /  tp  /  ckl , ck2 , ck3 , fOl , f02 , f03 , rfk( 5) , rcx(5) , alpha( 5) , 

&  c(5) , cone (5) , c three (5) , first (5) , second(5> , third (5) , fourth(5) , 

&  sevnth(5) ,elghth(5) .ctwo(5) ,flfth(5) , slxth(5) , tmld(5) , thlgh(5) , 
(t  ae  ,  rho(5)  ,xm(5,5),q(5,5)  ,e(5,5),z(5,5) ,  key  therm ,  theta  (5)  ,ck(5) 
common  /  be  /  tf lame . tf, epsilon, sigma, tc , tlm.dt , tbld, he , rt ,bi , 

&  chighsc(20), lay rout , Ire start , Icool , tlmcool ,uoa , tw 
If  (  keyeqn  . ne ,  0  )  go  to  30 

c  calculate  melting  energy  and  new  temperature  If  finishes  melting 

qm(l)  -  qm(l)  +  dt  *  (  (  ck2  +  .25  *  (  ckl  -  ck3  )  )  * 

&  (  rad(l)  +  .5  *  dx(ll)  )  *  (  t(l-l,l)  -  t(l,l)  )  • 

&  (  ck2  -  .25  *  (  ckl  -  ck3  )  )  *  (rad(l)  -  ,5  *  dx(ll)  )  * 

&  (  t(l.l)  ■  t(l+l,l)  )  )  /  <  rad(l)  *  dx(ll))  +  dt  *  dx(ll)  *  ae 

If  (  qm(l)  .gt.  hni2(ll)  )  go  to  10 
If  (  qm(l)  .It.  0.0  )  go  to  20 
t(1.2)  -  t(l,l) 
return 

10  t(l.2)  -  t(t.l)  4  (  qra(l)  -  hra2(ll)  )  /  (rho<ll)  ★  c(ll)*dx(ll)) 
qm(l)  -  hm2(ll) 
return 

20  t(l,2)  -  t(l,l)  +  qm(i)  /  (  rho(ll)  *  c(ll)  *  dx(ll)  ) 
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qm(i)  -  0.0 
return 

30  t(l,2)  -  -.(i.l)  +  (  (  f02  +  .25  *  (  fOl  -  f03  )  )  * 

&  (  rad(i)  +  .5  *  dx(ll)  )  *  (  t(l-l.l)  -  t(i.l)  )  - 
&  (  f02  -  .25  *  (  fOl  -  f03  )  )  *  (  rad(i)  -  .5  *  dx(ll)  )  * 

&  (  t(l,l)  -  t(i+l,l)  )  )  /  rad(l)  +  rfk(ll)  *  ae  *  (  dx(ll)**2  ) 
return 
end 

subroutine  tintface 

character*!©  matl(5) , explid.weapid, tltle(8) 
common  /  a  /  t(1000 , 2) , tm(5) , qm(lOOO) ,hm(5) ,hml ( 5) ,hm2 (5) ,k , 

&  kk(5),ll,12,13,m,l,pt(5), ipt(5) , rad (1000) , d(5) .keyeqn, icrit . 

&  ctr , maxt , dx ( 5) . dam(5) .n, Ipt .weapld . explid , mat! , Inc , r , title . 

&  V, time , tO , j j j ( 16) 

common  /  tp  /  ckl , ck2 , ck3 , fOl . f02 , f03 , rfk( 5) , rcx( 5 ) . alpha( 5 ) , 

&  c(5) .cone (5) ,c three (5) ,fLrst(5) ,second(5) , chird( 5) , fourth(5) , 

&  sevnth ( 5) , clghth( 5) , ctwo(5) .fifth (5), sixth (5), tmid(5) . thlgh(5) , 
6  ae , rho (5).xm(5,5).q(5.5).e(5.5),z(5.5) .key therm. theta ( 5) , ck(5) 
common  /  be  /  cf lame , tf , eps lion , sigma . tc , tim, dt . tbld .he , rt . bl , 

&  thlghsc(20) , lay rout . ires tart. Icool . tlmcool .uoa , Cw 
f(x , a.b.c ,d, e)  -  (  a  +  b*x  +  c*x*x  +  d*x*x*x  +  e*(x**4)  )  /241.9 
g(x,a.b,c,d)  -  (  a  +  b*(  1.  -  exp((c-x)/d)  )  )  /  241.9 
zz(x,a.b,c)  -  (  a  +  b/(x*x)  +  c/(x**4)  )  /  241.9 
c  ***  calculate  thermal  and  heat  generation  properties  ************** 
tbfor  -  (  (  t(m,l)  +  t(m+l.l)  )  /  2.  )  *  l.S  •  460. 
if  (  tbfor  .  ge .  tmld(ll)  )  go  to  5 

ck4  -f( tbfor .cone (11) . first (11) .second(ll) .third (11) , fourth(ll)) 
go  to  20 

5  if  (  tbfor  . ge .  thlgh(ll)  )  go  to  10 
ck4  -  22  (  tbfor,  ctwo(ll).  fifth(ll),  slxth(ll)  ) 
go  to  20 

10  ck4  -  g(tbfor  .cthreedl)  ,se\-nth(ll)  .thighdl)  ,clghth(ll)) 

20  aebfor  —  0. 
do  30  ii-l,kk(ll) 

30  aebfor  -aebfor  +  xm(ll,ll)  *  q(ll,il)  *  z(ll,li)  * 

&  exp( -e(ll , ll)/(r*tbfor) ) 
aebefor  -  aebefor  *  rho(ll) 

tafter  -  (  (  t(m+l,l)  +  t(m+2,l)  )  /  2.  )  *  1.8  -  460. 
if  (  tafter  . ge .  tmid(12)  )  go  to  35 

ck0-f( tafter .cone (12) ,first(12) ,second(12) , third (12) , fourth ( 12 ) ) 
go  to  50 

35  If  (  tafter  . ge .  thlgh(12)  )  go  to  40 
ckO  -  zz  (  tafter,  ctwo(12),  flfch(12),  sixth(12)  ) 
go  to  50 

40  ckO  -  g( tafter ,cthree(12) , sevnth(12) , thlgh(12) , elghth(12) ) 

50  aeaftr  —  0. 
do  60  li  -  l,kk(12) 

60  aeaftr  -  aeaftr  +  xm(12,li)  *  q(l2,ii)  *  z(12,ii)  * 

&  exp( -e(12 , il)/(r*tafter)) 
aeaftr  -  aeaftr  *  rho(12; 

c  ***  branch  to  proper  equation  to  calculate  new  temperature  ******** 
if  (  t(m+l,l)  ,lt.  tm(12)  )  go  to  90 
if  (  qm(m+l)  .gt.  hm(12)  )  go  to  100 
c  ***  calculate  melting  energies  ************************************ 
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qm(m+l)  -  qm(m+l)  +  dt  * 

&  (ck4  *  (rad(m+l)  +  .5  *  dx(ll))  *  (t(m,l)  -  t(m+l,l))  /  dx(ll)  - 
&  ckO  *  (rad(m+l)  -  .5  *  dx(12))  *  -  t(ra+2,l))  /  dx(12) 

&  +  (  rad(m+l)  •  .25  *  dx(12)  )  *  dx(12)  *  .5  *  aeaftr  )  /  rad(m+l) 
if  (  qm(ra+l)  .gt,  hml(12)  )  go  to  70 
if  (  qm(ra+l)  .It.  0.0)  go  to  80 
C(ni+1,2)  -  C<m+l.l) 
return 

c  ■***  calculate  temperature  rise  from  excess  energy  above  melting  energy 
70  t(m+l,2)  —  t(m+l,l)  +  (  qm<m+l)  -  hral(12)  )  *  rcx(ll)  /  dt 
qm(m+l)  -  hml(12) 
return 

c  **★  calculate  temperature  drop  from  resolidlf icatlon  ****★**********•*: 
80  t(m+l,2)  —  t(m+l,l)  +  qm(m+l)  *  rcx(ll)  /  dt 
qm(ra+l )  -  0  . 
return 

90  if  (  qm(m+l)  .le.  0.0  )  go  to  100 

Q  calculate  melting  energy  ^^★★★★★★★★★★★ii^'*’^**'*'*!*'*'**^^*^*****'**^******* 

qm(ni+l)  -  qm(m+l)  +  dt  * 

&  (ck4  *  (rad(m+l)  +  .5  *  dx(ll))  *  (t(m.l)  -  t(m+l,l))  /  dx(ll)  - 
&  ckO  *  (rad(mfl)  -  .5  *  dx(12))  *  (t(m+l,l)  •  t(m+2,l))  /  dx(12) 

&  +  (  rad(m+l)  •  .25  *  dx(12)  )  *  dx(12)  *  .5  *  aeaftr  )  /  rad(m+l) 
if  (  qm(m+l)  .It.  0.0  )  go  to  80 
if  (  qm(m+l)  ge.  hml(12)  )  go  to  70 
t(m+l,2)  -  t<m+l,l) 
return 

c  ■*■**  calculate  temperature  rise  (no  melting  or  re-solidiflcatlon)  **** 
100  t(m+1.2)  -  t(m+l,l)  +  rcx(ll)  * 

&  (  ck4  *  (rad(m+l)  +  .5  ★  dx(ll))  ★  (t(m,l)  •  t(m+l,l))  /  dx(ll) 

&  -ckO  *  (t(m+l,l)  -  t(m+2.1))  *  (rad(m+l)  •  .5  *  dx(12))  /  dx(12) 

6,  -f  (  rad<m+l)  -  .25  *  dx(12)  )  ♦  aeaftr  *  .5  *  dx(12)  )  /  rad(m+l) 
return 
end 

subroutine  prntrslt 

character*10  matl(5) ,explid,weapid, tltle(8) 
dimension  b(16),pos(16) 

common  /a  /  t(1000,2) , tm(5) ,qm(1000) ,hm(5) ,hml(5) ,hm2(5) ,k, 

&  kk(5) ,ll,]2,13,m,i,pt(5),ipt(5) , rad (1000) , d(5) .keyeqn, icrit , 

&  ctr ,maxt ,dx(5) ,dam(5) ,n, Ipt ,weapld,explld,matl , inc , r , title , 

&  v.time,tO,JJJ(16) 

common  /  tp  /  ckl , ck2 , ck3 , fOl , f02 , f03 , rfk(5) , rcx(5) , alpha(5) , 

&  c(5)  .  cone (5) , c three (5) , first (5) , 8econd(5) , third (5) , fourth<5) , 

&  sevnch(5) ,eighth(5) ,ctwo(5) ,fifth(5) .sixth (5) ,tmld(5) ,thlgh(5) , 

&  ae , rho (5),xm(5,5),q(5,5),e(5,5),z(5,5) .key therm, theta (5) , ck(5) 
common  /  be  /  tf lame , tf , epsilon, sigma , tc , tim, dt , tbld, he , rt ,bl , 

&  thighs t( 20) , layrout , ire start , Icool , timcool ,uoa , tw 

c- . output  point  control  J(l)  thru  J(16) . 

if  (  ctr  .ne.  0  )  go  to  167 

do  839  m6-2.16 
JJJ(n>6)-Jjj(m6-l)  +  inc 
8  39  continue 

wrice(6.1)  (title(ilp),ilp-l,8) 

1  format(8alO) 


write(6,10)  explid.weaptd.rt 

write(6,18)  dt 

wrice(6 ,106) 

wrice(6 , 17) 

write (6 , 23) 

do  171  i9-l,k 

write (6. 19)  19,matl(i9) ,dam(19) .pt(19) .dx(19) ,rho(19) ,c(i9) . 

lck(19) ,alpha(i9) ,theta(i9) ,tm(i9) .hm(19) 

171  continue 
write (6 , 24) 

24  format(lh0. 70hconductivtty  equations  (temp,  coeff.  units  in  "f  , 
lck(i)  units  in  "k)) 
do  172  i7-l,k 

write(6,32)i7  ,cone(i7) , first (17) ,second(i7) , third (17) , fourth (17) 

172  continue 

32  format(lh0,3hck(.i2.5h)  -  (,fl0.4,2h  - .el3 . 5 .4h*t  +.el3 . 5 , 6h*t*t  - 
l,el3.5.8h*t*t*t  +.el3.5.16h*(t**4))  /  241.9) 

wrlte(6 ,26) 

26  format(///20h  boundary  conditions) 
write (6 , 34) 

write (6 , 35) tO , epsilon, he , tc 

35  format(lh  , 3x , f7 . 1 , 9x , f 5 . 2 , 8x, f8 . 5 , 6x, f 8 . 1) 

34  format( IhO , 12hlnitial  temp,5x,7hepsilon, 

15x,  12hconvec  coeff ,  5x ,  9hcrlt  tetnp/3x ,  7h(deg  k),19x, 
215h(cal/sec-cm2*k) ,7x,3h(k)) 
if  (  icool  .ne.  1  )  go  to  341 
write  (6,342) 

342  format(///lh  ,12hTime  Cooling, 12x,8huoa,13x,6htw, 

6i/lh  ,12hStarts  (sec) , 9x ,  15h(cal/cm2-sec-K) , 9x, 7h(deg  K) ) 
wrlte(6,343)  timcool ,uoa , tw 

343  formatdh  ,  3x ,  f5 . 1 , 16x,  f8 . 6 , 13x .  f5 . 1) 

341  wrlte(6,33) 

write(6,36)  v,tlme,dt 
if(k.eq.l)  go  to  86 
do  85  ni2-l,13 
m3  -  m2  +  1 

85  write  ^  ,  3'’ )  [,.2. m3 

86  continue 

33  format(///  21h  interface  properties ,45x, 22hcomputatlon  parameters) 

36  format ( IhO , 3x , 8hlocatlon , 12x , 12hconductivlty , 30x , 4hv  »  ,fl0.6,2x, 
17htime  -  ,  f  7 , 2 , 2x  ,  5hdt  -  ,fl2.8/15h  (between  lyrs)  ,  12x ,  6hr.oef  f  . /) 

31  formatdh  ,  3x ,  2  ( 13) ,  18x ,  f 6 . 4) 
write (6 , 37 ) 
do  999  k9  -  l,k 

wrlte(6.  113)  k9,kk(k9) 

113  format(lh  , dayer  number  ' , il , 5x , 18hnumber  ind.  comps ., 19x , 13//) 
wrlte(6,107)  (e (h9 ,k5) ,k5-l ,kk(k9) ) 

107  format(lh0 . 30hactivaticn  energy  (cal/mole)  , 5 (lx , el5 . 5) ) 
write (6, 108)  (z(k9,k5) ,k5-l ,kk(k9) ) 

108  formatdh  ,  30hcollislon  number  (l./sec)  ,  5  ( lx ,  el5 . 5)  ) 

write (6, 109)  (xm(k9,k5) ,k5-l .kk(k9) ) 

109  format(lh  ,13hmass  fraction, 17x , 5(6x , f 10 . 5) ) 
writc(6,  110)  (q(k9,k5) ,k5-l,kk(k9)) 

110  format(lh  ,30hheat  of  reaction  (cal/gm)  , 5 (6x , f 10 . 5) ) 
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999  continue 

write(6 . 11) 
do  111  j7-l,16 
J6-jjj(j7) 

111  pos(J 7)-rad(l) -rad( j6) 

write(6,12)  (pos ( j 6)  , J 6-1 , 16) 

167  continue 

do  168  18-1,16 

168  b(18)-t(jl,l) 

write (6. 13)  tim.tf, (b(18) ,18-1.16) 
wrtce(9.1013)  n.tlm.maxt.raj 
wrtte(9.1014)  (t(l,2) ,1-l.n) 

1013  formatdh  ,  110 , 5x ,  f 7 . 3 , 5x .  110 . 5x,  110) 

1014  forroat(lh  ,8(lx,f9.5)) 
rewind  (9) 

11  format(62hlChe  following  table  Is  a  temperature  history  of  the  exp 
lloslve///) 

12  format<12h  depth  (cm. ) . 5x . 16 (2x. f5 . 2)//3x, lOhtime  free/2x . 13h(sec 
1.)  stream/9x, 5hflame/9x,  8htemp("k)/) 

13  formatClh  , f 7 . 2 , 2x, f6 . 1 , lx, 16 (lx, f 6 . 1) ) 

17  format(lh0, 127hlayer  material  thickness  pts  dx  density  s 
Ipcc.heat  therm,  con.  alpha  fourler  mod.  melt. temp,  heat 
2  of  fusion) 

23  formatdh  ,  18x .  5h(cm .  )  , lOx, 5h(cm.  )  . 2x,  7h(gm/cc) .  2s , 26h(cal/gm  "k)( 
Ical/sec  cm  "k) . llh(cm*cm/sec) ,18x,4h("k) ,9x,8h(cal/gm)/) 

10  format ( IhO , 2x , 29hexploslve  Identification, al0//2x, 

121hweapon  identification, 4x, al0//2x,9hradlus  -  , 

2f 10 . 6 , 4h  cm. ) 

18  format(lh0,2x, 16htlme  Increment  f7 . 5 , 2x, 3hsec) 

19  formatdh  ,  13 , 3x,al0 ,  lx,  £7 ,4,  lx,  f5 . 1 ,  lx,  f7 . 5 , 2x,  f7 .4 ,  lx,  2(3x,  £7 . 5, 
Ux  ),2(2x.f6.3,4x),2x,f7.2,7x,f7.3) 

37  format(///21h  explosive  properties) 

106  format(///20h  material  properties) 

return 

end 

subroutine  prntcrlt 
dimension  zx(16),b(16) 

character*10  matl(5) ,explld,weapld, tltle(8) 
common  /  a  /  t(  1000 , 2 ) ,  tm( 5)  ,qm(1000)  ,hm(5)  ,hml(5)  ,hni2  (5)  ,k , 

&  kk(5) .11,12. 13,m, l,pt(5) ,lpc(5) .raddOOO) ,d(5) .keyeqn, Icrlt, 

&  ctr ,maxt .dx(5) ,dam(5) ,n. Ipt . weapld.explld.matl , Inc, r , title , 

&  V, time , tO, J j j (16) 

common  /  tp  /  ckl , ck2 , ck3 , fCl , f02 , f03 , rfk(5) , rcx( 5) , alpha(5) , 

&  c(5) ,cone(5) ,cthree(5) ,flrst(5) , second(5) , thlrd(5) , fourth(5) , 

&  sevnth(5) ,elghth(5) ,ctwo(5) , f lfth(5) , slxth( 5) , tmld(5) , thlgh(5) , 

&  ae,rho(5) ,xm(5,5) ,q(5,5) ,e(5 , 5) . z(5 , 5) , keytherm , theta(5) , ck(5) 
common  /  be  /  tf lame , tf, epsilon. sigma, tc , tlm, dt , tbld, he , rt ,bl , 

(x  chlghst  (20)  ,  layrouc  ,  Ires  tart .  Icool .  tlmcool ,  uoa ,  tw 
do  10  jJ-1,16 

jl-JjJ(Jj) 

zx(JJ)-t(Jl,l) 

10  b(Jj)-t(jl.2) 

tlmb-tlm-dt 
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write (6 , 20)  timb , tf , (zx(i8) , 18-1 , 16) 
wrlte(6,20)  tim , tf . (b( i8) , t8-l . 16) 
write(6,30)  rad(k3) , tim, timb 
20  format(lh  , f7 . 2 , 2x, f6 . 1 , lx, 16(lx , f6 . 1) ) 

30  format(///65h  a  critical  temp  has  been  reached  within  the  explosiv 
le  at  radius-, f 5 . 2 , 18h  and  between  time-,f7 .2,4h  sec.lOh  and  time-, 
2f7.2,4h  sec) 
return 

end 
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APPENDIX  C 
PROGRAM  HTCOEFF 
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program  htcoeff 

dimension  tcool(5) ,  templ(5) ,  temp2(5),  uoa(5) 
character*!  answer 
character*15  fileS,  file6 

c  ***  layer  1  Is  the  steel,  layer  2  is  the  FM-26  (if  app)  ************** 
data  rhocpl  /  1.069  /,  ckl  /  0.109  /,  ck2  /  3.5E-04  / 
data  drl  /  O.A32  /.  dr2  /  0.406  / 
data  tcool  /  63.,  91.,  21.,  113.,  -118.  / 

write(*, 5) 

5  formate  Enter  name  of  input  data  file:  '\) 
read(* , ' (a) ' )  fileS 

write(*,6) 

6  formate  Enter  name  of  output  data  file:  '\) 
read(* , ' (a) ' )  filed 

open(5 , file-fileS) 

open(6 , file— filed , status-' new' ) 

write(*,8) 

8  formate  Is  this  case  for  a  calorimeter  with  FM-26?  (y  or  n)  '\) 
read(* , ' (a) ' )  answer 
write  (*,  7) 

7  formate  Enter  the  number  corresponding  to  the  coolant:  '/ 

&  '  1:  Water  fog,'/'  2:  Streaming  water,'/'  3:  AFFF,'/ 

&  '  4;  Halon  2402  or'/'  5:  Liquid  nitrogen') 

read  (*,*)  tz 
icount  -  0 

if  (  icount  .ne.  0  )  go  to  20 
read(5,*)  tirael,  (templ(i),  i-  1,5) 

10  icount  -  icount  +  1 

20  read(5,*, end-60)  time2,  (temp2(i),  i-  1,5) 

c  ***★  calculate  overall  heat  transfer  coefficient  ******************** 
do  40  i  -  1,5 

dtdt  -  (  terapl(i)  -  temp2(i)  )  /  (  timel  -  time2  ) 

q  -  rhocpl  *  dtdt 

if  (  answer  .eq.  'y'  )  go  to  30 

ts  -  q  /  (drl  *  ckl)  +  terap2(i) 

go  to  40 

30  tint  -  q  /  (dr2  *  ck2)  +  temp2(i) 
ts  -  q  /  (drl  *  ckl)  +  tint 
40  uoa(i)  -  q  /  (  ccool(lz)  -  ts  ) 

c  ***■  print  out  answers  *********************************************** 
write(6,45)  time2,  (  uoa(i)  , i-  1,5) 
writfe(*,45)  time2 ,  (  uoa(i)  , i-  1,5) 

45  formate  '  , 3x . f4 . 0 , 5 (9x, el2 . 6) ) 

c  ***  reset  variables  for  next  tinie  step  ****************************** 
do  50  J  -  1,5 
50  temp!(j)  -  cemp2(j) 
timel  -  time? 
go  to  20 
60  stop 
end 
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APPENDIX  D 

CALORIMETER  TEMPERATURE/TIME  PROFILES 
FOR  VARIOUS  COOLANTS 
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Figure  D-1.  Calorimeter  Temperature 
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Figure  D-2.  Calorimeter  Temperature/Time  Profile  (Uncoated,  Water  Coolant,  341  L/mln) 
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Figure  D-3.  Calorimeter  Temperature/Time  Profile  (Uncoated,  Water  Coolant,  170  L/min) . 
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Figure  D-4.  Calorimeter  Tempera ture/Tlme  Profile  (Uncoated,  Water  Coolant,  170  L/min). 
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Fig-re  D-5.  Calorimeter  Temperature/Time  Profile  (Uncoated,  AFFF  Coolant,  341  L/mln) . 


450 


Calorimeter  Temperature/Time  Profile  (Uncoated .  AFFF  Coolant,  341  L/mln) . 
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Figure  D-7.  Calorimeter  Temperature/Time  Profile  (Uncoated,  AFFF  Coolant,  170  L/min) . 
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gure  D-8.  Calorimeter  Temperature/Time  Profile  (Uncoated,  AFFF  Coolant,  170  L/min) 
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Calorimeter  Temperacure/Tlme  Profile  (Uncoated,  Halon  2402  Coolant, 


500 


OOOOOOOOOO 

loomoiooioolo 

'^'^COCOWCVJ'-t^ 


(o)  aanj,vHad?iai 


67 


Figure  D-IO.  Calorimeter  Temperature/Tine  Profile  (Uncoated,  Halon  2402  Coolant,  2.5  kg/s). 
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Figure  D-11.  Calorimeter  Temperature/Tirae  Profile  (Fii-26  Coating,  Water  Coolant,  341  L/min) 
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2.  Calorimeter  Temperature/Time  Profile  (FM-26  Coating,  Water  Coolant,  3A1  L/min) 
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(Th«  reverse  of  this  pago  Is  blank.} 


Figure  D-16.  Calorimeter  Temperature/Time  Profile  (FM  -26  Coating,  Liquid  Nitrogen  Coolant) 


